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Abstract 

We describe the investigation of spontaneous mass-generation and 
chiral symmetry breaking in supersymmetric QED3 using numerical so- 
lutions of the Dyson- Schwinger equation together with the CJT effective 
action and supersymmetric Ward identities. We find that, within the 
quenched, bare vertex approximation, the chirally symmetric solution is 
favoured. 

1 Introduction 

Any realistic attempt to unify gravity with the other fundamental forces of 
nature must almost certainly be supersymmetric. Whilst perturbative re- 
sults are often easier to find in supersymmetric field theories than in ordinary 
gauge theories, especially with the use of superfields, the development of non- 
perturbative tools for supersymmetric theories has not yet reached maturity. 
In this paper we apply a nonperturbative tool which has proved to be effec- 
tive in .the study of strong interactions, namely Dyson-Schwinger equations 
(DSEs)tH, to one of the simplest supersymmetric gauge theories. 

We choose as our gauge theory supersymmetric quantum electrodynamics 
in three space-time dimensions (SQED3). In its nonsupersymmetric form, 
three dimensional quantum electrodynamics (QED3) with four component 
spinors □ exhibits the interesting nonperturbative phenomena of confinement 
and, at least for small numbers of flavours, chiral symmetry breaking. Much of 
what js known about the behaviour of QED3 has been learnt from numerical 
DSE El or lattice □ calculations. Considering the success of DSE techniques in 
nonsupersymmetric theories, there is a surprising scarcity of applications of 
the method to supersymmetric theories in the literature, and an almost total 
absence of numerical DSE calculations in particular. 

The version of SQED3 whiclxwe consider is the four-component fermion 
version first proposed by PisarskiB, who approached the model by dimensional 
reduction from SQED4. Herein, we develop SQED3 using the Wess Zumino 
construction a. Our construction produces supersymmetry multiplets which 
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differ slightly from those of Pisarski because of an extra degree of freedom in 
the charge conjugation matrix. Pisarski's approach was an analytic one based 
on a 1/iVflavor expansion. His analysis indicates the existence of a dynamical 
mass generating solution in the large A% avo r limit. 

Koopmans and SteringaB have analysed SQED3 with 2-component fermions 
using DSEs and claim the critical number iV c = 3.24 of flavours below which 
spontaneous mass generation occurs. This result should be treated with some 
caution however as it ignores the possibility of spontaneous parity violation 
and the generation of a photon mass via a Chern-Simons term when the num- 
ber of two-component fermions is odd. Furthermore, there is no indication 
whether spontaneous mass generation corresponds to the stable vacuum of the 
theory. 

The above treatments use a component formalism to study SQED3. By 
contrast, Clark and LoveB develop DSEs for SQED4 using a superfield formal- 
ism. While more elegant as a mathematical formalism, the superfield approach 
has the disadvantage that each DSE contains an infinite number of terms. This 
is dealt with by truncating from the chiral multiplet self energy diagrams con- 
taining seagull and higher order n-point vertices. They find that the effective 
mass contains a prefactor which vanishes in Feynman gauge and conclude that 
there can be no spontaneous mass generation in SQED4. 

The work of Clark and Love has been criticised by Kaiser and Selipskyi 
on two grounds. Firstly they argue that the truncation of seagull diagrams is 
too severe as it ignores contributions even at the one-loop level. Secondly they 
point out that infinities arising from infrared divergences which plague the su- 
perfield formalism can counter the vanishing prefactor and allow spontaneous 
mass generation. These criticisms highlight some of the dangers of attempting 
to extract phenomenological consequences of supersymmetric DSEs by work- 
ing solely .-with the superfield formalism. In fact, analyses in the literature 
of SQCD 113, have generally found the component formalism to be the most 
efficient way to proceed. 

In this paper we present a first detailed numerical analysis of the chiral 
multiplet DSE for a supersymmetric gauge theory. We work with the com- 
ponent formalism rather than the superfield formalism to avoid the problems 
encountered by Clark and Love. By choosing to work in three dimensions, 
and with the component formalism, we find that all integrals encountered are 
infrared and ultraviolet convergent. 

Section || describes the four-component fermion Clifford algebra of QED3 
and gives the supersymmetry multiplets. In section ^ we present the SQED3 
lagrangian and chiral multiplet propagator DSEs. We explain how supersym- 
metric Ward Identities can be used to simplify what would otherwise be an 
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arduous problem and present our numerical solutions for quenched SQED3 
in the bare vertex approximation. We present both chirally symmetric and 
asymmetric solutions and use the Cornwall- Jackiw-Tomboulis (CJT) effective 
potential to determine that the massless solution is dynamically favoured. 



2 The algebra of SQED3 

In 2+1 dimensions there are two inequivalent, irreducible representations of 
the Clifford algebra, given iiLterms of 2 x 2 matrices. These two representa- 
tions differ by a minus signlHI, and have the undesirable property that either 
one leads to a version of QED3 which is parity non-invariant. To circumvent 
this property, it is common to consider a four-component version of QED3 
incorporating Dirac matrices which are a direct sum of the two inequivalent 
representations □. The Dirac matrix algebra we employ here is constructed as 
follows 113: 

The 4x4 matrices 7 M satisfy {7^, 7^} = 2rj^, rj^ = diag(l, -1, -1) where 
/1 takes the values 0,1 and 2. We take the complete set of 16 matrices 

{1a\ = {^,74,75,745,7^7^4,7^7^45}' 
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m = - = 1 j Q j , 75 = 7 = y u q )■'.!-. = ' r..-,. 

7^4 = «7m74, 7m5 = HfJ.15, 7M5 = 7m 745 • 

The matrices 1,74,75,745 are the Pauli matrices in block form and as such 
generate a U(2) algebra. The parity and charge conjugation rules for four- 
component Dirac spinors are given by 

Vip{x)p- 1 = n^(x°,~x\x 2 ), V^{x)V- x = $(x°, -x\ x 2 )^ 1 , 

CipC- 1 = Ctp T , C$C~ l = -^ T c _1 , (2.1) 

where 

II = 7i 4 e^ P74S , C = 72 e^ c745 , (2.2) 

and (0 < (j)p,(j>c < 2tt). A Majorana spinor is one for which ip = Cip T . The 
arbitrary phases <pc and 4>p are important for classifying the bound states in 
QED3tl. 
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We have 



n-M 74 )ii = r p 74 



(2.3) 



where 



_ / - cos 20 P - sin 20p \ _ / - cos 20 c sin 2<f> c \ ,„ 
p \-sm2^)p cos20 P )' c \ sin 20 c cos20 c /' { ' 

The first step in the extension of QED3 to supersymmetry is the construc- 
tion of chiral multiplets. A reasonable tentative suggestion is 

5a = —iCijj 
5b = (j 5 ip 

H = (/ + «75<?) + i $>{a + il 5 b)( (2.5) 
Sf = 

5g = <75 

where a, b, f and g are real and £ and if> are Majorana. If the superalgebra is 
to hold then the commutator of two of these transformations must obey 

[6 l ,h]X = 2C2-fC l dpX, (2.6) 
where X is any component of the multiplet. For the case of a we have 

[5 1 ,5 2 }a = frf^a - i&TlsVW -(!«—> 2). 
Now, since £ is Majorana 

Ci7 M c 2 = -Grci, (2.7) 



as required, but using equations (|2.3|) and (2.4) 

-<i757 p C2 = C2(75Cos20 c + 74sin20 c -)7' 1 Ci ) (2.8) 

which does not cancel C2757 M Ci- A similar situation arises if a is replaced by 
any other member of the multiplet. To stop this 'blowing out' of terms we could 
simply set <pc = 0. This would be unfortunate though as the angular freedom 
<t>c in the matrix C would be lost. As stated earlier, this angle is important 
for classifying the bound states in QED3 so we would like to preserve it if 
possible to see what effects it has (if any) in the supersymmetric theory. To 
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this end we define the rotated Dirac matrices. This can be done by making 
the substitution 

74 \ Hp \ = (cosfc -sm0c\ (lA =M (lA {2 g) 

15 J \1W J V Sm< PC COS( PC J \75 J \l5 J 

in the Clifford algebra. (Note that — ijpjw = 745 so the matrices /, 7p, jw, 745 
again generate an SU(2) algebra.) Then 

c -0:)^(-l)' < 2 ' io » 

since 

MR C M- I =(~*°y (2.11) 

The rotated matrices can be used to define a supersymmetry transformation 
consistent with Eq. (2.6). In order to do this, all terms except Ci7 A 'C2 generated 
by 6162X must be symmetric under interchange of Ci, C2 ; i-e. 

blwl^Ci = Ci7iv7 M C 2 . (2.12) 
In this sense ~/w is well-behaved but jp and 745 are problem matrices because 

C 2 (7P,74 5 )7 M Ci = ~Ci(7P, 745)7 M C2- (2.13) 
Making the substitution ( |2.9| ) gives the 2+1 dimensional chiral multiplet 

5a = 

Sb = (jwip 

H = {f + Hwg) + il^d^a + i lw b)C (2.14) 
Sf = (JtP 

which is the analogue of the standard chiral multiplet in 3 + 1 dimensions. 

The relative difference between the arbitrary phases <pp and <pc is fixed by 
the imposition of supersymmetry. Indeed, from equations (2.3) and (2.£) we 
have that 

IT 1 ( lP \li = MRpM- 1 ( 1P 



where 



MR P M-i = [ cos2 Jii P tt c \ " Sin 9 r?^ r ) • <*•!■•' 

1 -sin2(0p + 4>c) cos2(0p + c ) ' 



If the form of the chiral multiplet transformation Eq. (2.14) is to be maintained 
under parity transformations, the off-diagonal terms in this matrix must be set 
to zero. <fip must therefore be set to one of 



4>ck - 4>c 



(2.16) 



and we choose the former. The Clifford algebra to be used is now 



With the Clifford algebra and chiral multiplets established we now look 
for a general multiplet. This is very similar to the standard general multiplet 
in 3+1 dimensions except for an extra scalar held K, which makes up the 
bosonic degree of freedom lost when the vector field is taken from 3+1 to 2+1 
dimensions. Our general multiplet V is defined by the following fields and 
transformations: 



5C = ("fwX 

S X = (M + iywN)( + i^(A„ + i-y w d„C)( - lP K( 
SM = C(^X + «A) 

5N =iC7w(ftC+*A) (2] 

5 An = C7m_ a - K^nX 
5K = — i^jpX 

S\ = I( 7 - 7 m _ ^Y)d„A„( + i lw D( + i lP $K(, 
SD = iClw^X. 

3 Solving the DSEs of SQED3 

Repeating the method of Wess and Zuminoi yields the SQED3 lagrangian 



1A = U,7p, 7^,745,7^,7/^,7/^, 7^5} 
C = 72 e #c745 , n = 7l7P . 



(2.17) 



L=\f\ 2 + \g\ 2 + \d^ + \d^-^ W 



m{a*f + af* + b*g + bg* + iipip) 



ieA»{a* d^a + b* d^b + ^ 7m V) + eK^jpifj 



(3.1) 



e[X(a* + ijwb*)ip — ip(a + i-fwb)X 



+ieD{a*b - ab*) - e 2 (K 2 - A^)(\a\ 2 + |6| 2 ) 



\F^ V F L 
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|A @\ + \d»Kd^K + \D 2 
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which becomes that found by PisarskiB by dimensional reduction of SQED4 
when <f>c is set to zero and the scalar fields are trivially redefined. 

The chiral limit is defined by taking m —* 0. In this limit the bare la- 
grangian is invariant with respect to a global U (2) symmetry generated by /, 

74, 75 and 745. In the non-supersymmetric chiral theory, spontaneous mass 
generation leads to a nonperturbative breaking of the chiral generators 74 and 

75. Here we shall explore the possibility of chiral symmetry breaking in SQED3 
by considering the chiral multiplet propagator DSEs. We shall use Witten's 
result that supersymmetry remains unbrokenEjt 
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Figure 1: The Dyson-Schwinger equation for the electron propagator in SQED3. 



The electron propagator DSE in SQED3 is more complicated than that in 
the non-supersymmetric version since the electron interacts not only with the 
photon, but also with the photino, the K and its own superpartners a and b. 
The DSE for the electron is illustrated in Fig.|l], where a solid line indicates an 
electron propagator, dashed line an a or & propagator, a wavy line a photon 
propagator, an alternating dash-and-dot line a K propagator and a wave with 
a solid line through it a photino propagator. 

The DSE for the scalar partners is shown in Fig.|[ It would appear prima 
facie that we have two coupled DSEs to solve, the second considerably more 
complicated than the first. We are saved from ihis effort however, by the 
existence of supersymmetric Ward Identities (WI)EJ. The supersymmetric WI 
relating the scalar two point functions to the electron propagator is 



(^} = i{a*f) - i Aa*a) = i(b*g) - i ^{b*b). 
Substituting in the fermion propagator Ansatz 



(3.2) 



5(p) = (M) = 



i>A( v 2 ) + B{p 2 ) 



(3.3) 
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Figure 2: The Dyson-Schwinger equation for the a and b propagator in SQED3. 



gives the scalar propagator 

D{p 2 ) = (a*a) = (b*b) 
and the a to f amplitude 

(a*f) = (b*g) = 



A(p 2 



p 2 A 2 (p 2 ) - B 2 {p 2 Y 



_ B{p 2 ) 

p 2 A 2 {p 2 ) - B 2 {p 2 ) A(p 2 ) 



B(p 2 ) 



D( P 2 



(3.4) 



(3.5) 



If nonperturbative chiral symmetry breaking is allowed for, the two point func- 
tion (a* f) is potentially non-zero in the chiral limit even though perturbation 
theory predicts it should be identically zero. 

Clearly it suffices to solve only the fermion DSE to determine the scalar 
functions A(p 2 ) and B(p 2 ). To do the calculations we perform a Wick rotation 
to Euclidean space. Our Euclidean conventions are the same as those set out 
in Appendix B of reference Il3. Spacelike momenta satisfy p 2 > 0, and the 
fermion and boson propagators take the form 



S(p) 
D(p) 



17 -pA(p 2 ) +B(p 2 )'' 
A(p 2 

p 2 A{p 2 ) + B 2 {p 2 Y 



(3.6) 
(3.7) 
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Figure 3: The chirally asymmetric solution for the propagator function B(p 2 ) in Feynman 
gauge, £ = 1. If chiral symmetry is not broken, B is identically zero. 



respectively. We use the chirally quenched approximation^ which is known to 
be compliant with U(l) gauge invariance, that is, we take the propagators of 
the photon and its superpartners to be bare. The class of covariant gauges, for 
which the photon propagator is given by 



fc 4 



(3.8) 



will be considered. The corresponding K and photino propagators are 

D K (k) = ±, (3.9) 
1 



D x (k) = - 

ij ■ k 

respectively. We also employ the bare vertex or rainbow approximation 



(3.10) 



Fp(g,p) = -ierf P , 



(3.11) 
(3.12) 
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A(p) 3 




Figure 4: The chirally asymmetric (solid curve) and symmetric (dashed curve) solutions for 
the propagator function A(p 2 ) in Feynman gauge, f = 1. 



(3.13) 
(3.14) 



for the Euclidean electron-photon, electron- K and electron-a-photino vertices 
respectively. The Euclidean space DSE with these approximations reduces 
after angular integration to the following coupled integral equations for B{p 2 ) 
and A(jp 2 ): 



B(p 2 ) = (£ + 3) 

„2 



dq 



qB(q 2 



4m 2 p- jo 



dq 



dq 



q 2 A 2 {q 2 )+B{q 2 
qA{q 2 ) 



In 



P 



(3.15) 



q 2 A 2 {q 2 ) + B(q 2 ) 



p2 _|_ g2 



2p 



hi 



qA(q 2 ) 



2ir 2 pJ H q 2 A 2 {q 2 ) + B(q 2 ) 



In 



p + q 



p-q 



(3.16) 



We solve equations ( 3. 15 ) and (3.16) numerically using the standard itera- 
tive procedure introduced by Applequist et alQ. The functions A{p 2 ) and B{p 2 ) 
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Figure 5: The chiral condensate of quenched SQED3 (solid curve) and QED3 (short dashes) 
in the bare vertex approximation. Also plotted is the chiral condensate of quenched QED3 
using the minimal Ball-Chiu Ansatz for the fermion-photon vertex (dashed-dot curve) and 
the chiral condensate of quenched SQED3 using the partial approximation to the Ball-Chiu 
Ansatz described in the text (long dashes). 



are denned on a non-uniform grid of fifty-one points concentrated at small mo- 
menta where the function varies more rapidly. The integrand is interpolated 
using a cubic spline using an ultraviolet cut-off of p = lOOOe 2 . 

We show in figures ^| and |^ both the chirally symmetric and asymmetric 
solutions to the massless (to = 0), Feynman gauge (£ = 1) electron DSE of 
Fig.|l|. The results are plotted in units with e 2 = 1. It is apparent in figure 
|3| that mass generation is remarkably suppressed and both graphs descend 
steeply to the values assumed when there is no dressing. We have tested the 
convergence of our numerical iteration procedure by varying the initial guess 
for A and B and the ultra-violet cutoff. These changes had no significant effect 
on the solution obtained. 

The DSE is the first functional derivatijie with the respect to the electron 
propagator of the CJT effective potential and its solutions are therefore 
stationary points. To ascertain if the chirally asymmetric solution to the DSE 
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is dynamically favoured we make use of the CJTL effective potential which, 
when evaluated at a stationary point, is given by 113 

V[S,D] = J&L{ t r ln[l - E s (p)S(p)} + ±tr[E s (p)S(p)]} 



(2tt) 3 
d 3 p 



(2n) 



{tr ln[l - S D (p)£>(p)] + fa\E D (p)D(p)]} ,(3.17) 



where E s (p) = S(p)- 1 - S^fo) and T, D (p) = Dip)- 1 - D^ le (p). This 
formula neglects the contribution from the photon and its superpartners since 
it is calculat ed in the quenched approximation. Note that there is another 
term in Eq. ( 3.17| ), given by 

d 3 p d 3 q 



{2 ^(2^ [Dfifliq)+DKiq)]D(pl (3 - 18) 
corresponding to the vacuum graphs giving rise to the tadpole contributions 



in the boson SDE. As we have written it, Eq. (3.17) only gives half of these 
particular diagrams. However_in the quenched approximation Eq. (3.1S) is 
zero by dimensional reductionEO. 
In Table g we list the difference 

v4s,n ] -v s[ s,D ] = 2j^{J p2A ^ 2 



p 2 4(p 2 ) 2 + 5a(p 2 ) 
P 2 A A (p 2 ) 1 



P 2 A A ( P 2)2 + B A (p2)2 ^ s( p2) 



(3.19) 



between the effective potentials calculated for the chirally asymmetric and 
symmetric solutions for a range of values of the gauge parameter £. Our cal- 
culation confirms that the asymmetric solution is dynamically favoured in the 
nonsupersymmetric case. For SQED3, however, we find that chiral symmetry 
is not spontaneously broken, at least in the quenched, rainbow approximation. 

It is clear from Table [l] that the calculated effective potential is dependent 
on the gauge parameter £. This is an artifact of the rainbow approximation 
which, for dressed electron propagators, violates the U(l) Ward-Takahashi 
identity (WTI). Another useful indicator of gauge symmetry breaking is the 
U(l) invariant chiral condensate 

W>=trS(* = 0) = -^ dp ^^^ , (3.20) 

calculated in the asymmetric phase. In Figure ^| we plot the chiral condensate 
for quenched rainbow SQED3 and QED3. While the calculated condensate for 
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Table 1: The difference V A [S, D] - V s [S, D] between the CJT effective potential of the chirally 
asymmetric solution and the chirally symmetric solution for QED3 and SQED3 at different 
values of the gauge parameter £. 





SQED3 


QED3 





4.62 x 10~ 3 


-1.32 x 10" 5 


0.5 


4.29 x 10" 3 


-6.02 x 10~ 6 


1 


4.07 x 10" 3 


-3.44 x 10" 6 



SQED3 is surprisingly insensitive to the choice of gauge, the condensate for 
QED3 is, as expected, strongly gauge dependent. 

For QED3 the invariance of the chiral condensate can be considerably 
improved by replacing the bare vertex with the minimal Ball-Chiu vertex 
AnsatzEa, which is specifically designed to respect the C7(l) WTI and be free 
of kinematic singularities For comparison, the QED3 condensate obtained 
in this way in reference Ej is plotted in Figure ^. Here we attempt a similar 
substitution for SQED3. Note that the photon's supersymmetric partners K 
and A are completely invariant under a gauge transformation and their ver- 
tices are not directly constrained by the WTI. Compliance with the WTI can 
therefore be achieved by replacing the bare photon-fermion vertex with the 
minimal Ball-Chiu Ansatz whilst the remaining vertices are kept bare. This 
method incurs the penalty of breaking supersymmetry. The resulting chiral 
condensate is plotted in Figure [|. Suprisingly the variation of the condensate 
with respect to the gauge parameter was found to be an order of magnitude 
greater than in the bare case. We attribute this to the violation of super- 
symmetry and conclude that any attempt to improve the vertex must remain 
supersymmetric. 

4 Conclusions 

We have used the method of Wess and Zuminoi to obtain a lagrangian for 
SQED3 with four-component fermions. The lagrangian derived in this way 
agrees with that obtained by dimensional reduction □ except that our form 
retains an angular degree of freedom associatecLwith the definition of parity 
and charge conjugation in (2 + 1) dimensions E3. This is effected by taking 
suitable linear combinations of the Dirac matrices 74 and 75. It was found 
that supersymmetry constrains the angles <pp and 4>c associated with parity 
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and charge conjugation transformations in QED3EI to be related by <j)p = —4>c 
or 7r — cf>c ■ 

The DSE of SQED3 was solved numerically in quenched, rainbow approx- 
imation after using supersymmetric Wis to relate scalar and fermion propa- 
gators within the chiral multiplet. Application of the CJT effective potential 
revealed that the chirally symmetric solution is dynamically preferred, in con- 
trast to the situation in QED3. We conclude that chiral symmetry is not 
broken in SQED3 in the quenched, rainbow approximation. 

Although the rainbow approximation violates the U(l) gauge symmetry, 
we find the chiral condensate to be relatively insensitive to the choice of gauge 
fixing parameter. We believe this to be fortuitous however, since replacing 
the bare vertex with an Ansatz which respects the U(l) WTI but not the 
supersymmetric Wis produces a chiral condensate which is significantly gauge 
dependent. Work on developing a vertex Ansatz which respects both the U(l) 
gauge and supersymmetric Wis is currently in progress. 
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